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Abstract 

Soft and virtual loop photonic contributions to the second order next-to-leading 
QED radiative corrections to Bhabha scattering are calculated with help of the 
renormalization group approach. The results are in agreement with earlier calcula- 
tions, where other methods were used. Scale dependence and the present theoretical 
accuracy of Bhabha cross section description are discussed. 
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1 Introduction 



Precision theoretical predictions for the differential cross section of Bhabha scattering are 
of ultimate importance for all experiments at electron-positron colliders. They are required 
for normalization purposes including luminosity determination, for several searches of new 
physics, and as a background contribution to many other processes studied at the colliders. 
To provide the required accuracy of the predictions, one should take into account radiative 
corrections in the first and higher order of perturbative QED. Certain effects of strong 
and weak interactions should be included as well. 



It is natural to expand the QED part of radiative corrections to Bhabha scattering into 
a series in the fine structure constant a and in powers of the so-called large logarithm 
L = ln(M^/mg), where M is a large energy scale related to the beam energy, M ^ rrie- 
The terms enhanced by the large logs are known to give the bulk of the result in the 
kinematical regions of interest for experiments. 
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In this paper we present the derivation of the second order next-to-leading virtual and 
soft photonic radiative corrections to the cross section of Bhabha scattering. To get the 
corrections we use the renormalization group techniques borrowed from QCD. The result 
is found to be in agreement with earlier calculations [1,2]. The advantage of our approach 
is its universality: in the same way one can get all the other remaining contributions 
to the radiative corrections in O {o?L) . For the case of small angle Bhabha scattering 
the complete result for the O {a^L) is known [3] , but for the general case of large angle 
scattering terms of this order are still not systematized, while there is a number of results 
for particular contributions scattered in the literature. 

The paper is organized as follows. The notation is introduced in the next section. The 

application of the structure function approach is given in Sect. 3. The numerical results 
arc presented in several plots followed by a short discussion and conclusions. In sake of 
completeness, some useful formulae are listed in the Appendix. 



2 PreliminEiries 

Let us represent the differential Bhabha cross section as a series in a: 

da = da^"-"" + d(7(^) + d(7(2) + O {ot') , (1) 

where d(T^°™ represents the Born level cross section, and dcr^^'^-* are the pure QED con- 
tributions of the first and second order corrections. 

The first order contribution is usually decomposed into three parts: 



da^^) = d(7^ + d(7^ + d(7«, (2) 

where superscripts "V", "S", and "H" are used to denote the virtual, soft, and hard 
photonic corrections, respectively. The small parameter A (A ^ 1) subdivides the kine- 
matical domain of real photon emission into soft and hard parts with photon energy below 
and above A ■ i?bcam, where fi'beam is the beam energy. The one-loop contributions are well 
known, see i.e. Refs. [4,1,5]. 

In the second order we construct a similar decomposition: 



where the superscripts have the same meaning as in Eq. (2), so that for instance "VH" 
denotes the contribution due to emission of one hard photon accompanied by the effect 
of a single virtual loop. 
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In this paper we don't take into account the so-called pair contributions, related to emis- 
sion of real or virtual pairs {e~^e~ , t:^t:~ etc.). Their numerical contribution to the observed 
cross section is typically small compared to the photonic correction, see Ref. [6,7,8]. 

The second order contribution dcr^^^ is decomposed into series in the powers of the large 
logarithm. There in particular contributions we meet terms of the orders O (q;^L^'^'^'^'°). 
The terms with the fourth and the third powers of L will cancel out in the sum of virtual 
and soft photon contributions. 

Using the factorization properties of soft photon radiation, we can write immediately the 
Soft-Soft (SS) and Soft- Virtual (SV) contributions [9]: 

where 5^'^ — dcr^'^/dcr^"'^'^ (see Appendix for explicit formulae). Note that for the contri- 
bution of double soft photon emission, do"^^, in the formula above we apply an upper limit 
on the energy of each of the photons independently. The Virtual- Virtual (VV) contribu- 
tion can't be received in such a simple manner. Below we will show how to reconstruct 
the logarithmically enhanced part of it using the renormalization group techniques or in 
other words the electron structure function approach. 



3 Structure Function Approach 

The structure (fragmentation) function approach, widely used in QCD can be applied 
to QED problems [10,11,12,13,14]. With help of it we can analytically find the most 
important contributions reinforced by the large logarithm L, since they can be treated as 
electron mass singularities. 

We are going to drop the pair contributions, so we need here the pure photonic part of 

the non-singlet structure (fragmentation) functions for the initial (final) state corrections. 
These functions describe the probability to find a massless (massive) electron with energy 
fraction z in the given massive (massless) electron. In our case with the next-to-leading 
accuracy we have 

+ (^)' (^^'^ie' ® Pil\z) + LPi°) d,{z, /.o, me) (5) 
where the superscripts "str" and "frg" are used to mark the structure and fragmentation 
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functions, respectively. The difference between tfie functions appear only due to the dif- 
ference in the next-to-order splitting functions p(^'T). The modified minimal subtraction 
scheme MS is used. Wc have chosen the factorizations scale equal to M, and the renormal- 
ization scale will be taken equal to mg. More details on application of the approach to 
calculation of second order next-to- leading QED corrections can be found in Refs. [11,14]. 

The function di corresponds to next-to-leading order initial condition in the MS scheme 
for iterative solutions of the QED DGLAP equation. It reads 



di{z,iJ,o,me) 



^fln4-21n(l-;.)-l 
1 — z \ mi 



(6) 



where the usual plus prescription is exploited to regularize the singularity at 2; — > 1. The 
relevant lowest order splitting function is 



1 + z 



2n 



1-Z 



(7) 



The next-to-leading order functions P^l''^^ is provided by the set of Feynman diagrams 
with pure photonic corrections. 



pa.7)frg(^) = ^(1 -z)(^^- 3C(2) + 6C(3)) + l-l^(^21n^ln(l - ^) - 21n'^ 



2Li2(l-^) ) + -{l + z)\n^z + 2z\nz-3z + 2, 

1 + z^ 



(8) 



p(l''r>'^{z)^S{l 3C(2) + 6C(3)) + -±^(-2\nz\n{l - z) +\n'' z 

+ 2Li2(l-^)) - ^(l + ^)ln2^ + 21n^ + 3-2^. (9) 

Note that in these splitting functions the parts proportional to the delta-function are 
equal, and the remaining functions are not singular at 2; — > 1. Therefore, since at the 
present moment we are interested only in the soft and virtual part of the corrections, we 
can forget about the difference and drop the superscripts distinguishing the structure and 
fragmentation functions. 

The master formula describing the radiatively corrected Bhabha cross section in the struc- 
ture function approach reads [5]: 



da^J dzi J dz2Vf:{zi)Vf^{z2) (da^°™(^i, z^) + da^'\zi, Z2) + O 
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X /f /f (10) 

yi 2/2 

where da^^^ is the O (a) correction to the massless Bhabha scattering, calculated using 
the MS scheme to subtract the lepton mass singularities. Energy fractions of incoming 
partons are Zi^2, and Y12 are the energy fractions of the outcoming electron and positron. 

Here we are interested in the contributions due to virtual and soft photons, so all the 
four integrals will have the same lower limit being equal to 1 — A. First we can perform 
convolution of the four structure functions entering Eq. (10) with each other: 



D^\z) = {6{1 -z) + ^d,{z) + ^LP('\z)) 

^5{l-z)+4^(d,{z) + LP^'Hz)^ +4l^-^y(2L'{P^'\z)r 

+ 4Ldi{z)®P^''\z)+LP^^\z)^ +0(a^L'',a^), (11) 

with a short notation for multiple convolution operation: (^4)®" = A ® A ® A ■ ■ ■ A. If 

n 

2; = 1 — A and A <C 1, this function gives the probability to find such a situation where 
one looses in total due to photon emission Ai^bcam from the total energy of the process 
under consideration (the center-of-mass reference frame is used throughout the paper). 

Let us fix now the factorization scale M — y/s and define 



Later on we will consider another choice of the scale. 

Convolution of the function found above with the Born part of the kernel cross section 
gives us the corresponding part to the cross section (with the upper limit on the energy 
lost): 



J P®^(^)da^°™(^)d;2 = da^^^jl + ^ 4L/21n A + ^"j + C 

l-A 



SLl (P(0))^' + 16L,(P0 cii)A + 4L,Pi^)] + O (a^L^ a') , (13) 



where we used subscript A to specify the the so-called A-part of the corresponding func- 
tion (see i.e. Refs. [12,13]). 
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Convolution with the da^^^ is more comphcated since the latter is a non-trivial function of 
z. But we restricted ourselves to consider only the terms reinforced by the large logarithm. 
Therefore we need to compute only the following part: 



1 

a 



A^Ls / d;2[p(°) ®da«](^). (14) 

ZTT 



1-A 



Using the techniques of dealing with the singular functions regularized by introduction of 
the A and parts [12,13,15] we cast it into the following form: 



'T. i d./|L.P<»'g)d.<..W^4i^d.i'>..Pr.4i^...,.a.). (15, 

1-A 

where 



± 

J ZTT 



1-A 
1 



5(a«)= / da^^\z)2 \nll-^—^)dz. (16) 



1-A 



Let us note that the last non-trivial term the function under the integral above is not 
singular for z — > 1. So the virtual loop contribution to the function da-^^^ doesn't contribute 
to the result of the integral. And we can put the upper limit of the integral over z to be 
equal to 1 — Ai, Ai <S A ^ 1. Therefore we need to consider only the soft photon part 
of the function: 



da(''^\z) 
dz 



= d(7 



I 47r2 dz ^\k\2 + X2 



2n 



AA I A t , Pe 



Pe 



-1 



v^+k pek p'Ji p'+k^ 



a 
2^ 



di(^)+L,pW(;2) 



(17) 



where 1 — z = |A;|/£'beam; A is a fictitious photon mass, A <^ 171^; Pe (Pe+) is momentum 
of incoming electron (positron), p'^, p'^+, and k are the momenta of the outgoing electron, 
positron, and the photon. Using the standard techniques of calculations of soft photon 
contributions we get 



6{a 



1-Ai 



I 

1-A 



dz 



dz 



2 In ( 1 - 



1-z 



6 



a f az 



27r 



It 



l-Ai 



27r 



-4C(2)ln 



81n(l +41n^ - 4ln^-^ 



+ 8 C(3)-C(2)lnA 



l + c 



(18) 



where c is the cosine of the electron scattering angle, c = cospeP'e- 

Summing up the contributions in Eq. (15) and then summing up Eqs. (4,13,15) we receive 
the leading and next-to-leading second order contributions to quasi-elastic Bhabha cross 
section, where the total energy loss (due to soft photon emission) is limited by AE'bcam- It 
is useful to describe also the case when the energies of the soft photons (if they are two) 
are limited independently. The transition between the two cases was derived in Ref. [9]. 
Applying it we get the final result: 



da^v + da^^ +da' 



sv 



,ss 



a \ 
2^) 



321n^A + 481nA+ 18 



64 ( In 



X 



1 In^A 



+ 16(2Li2(l-x) -2Li2 (x) + 31n(^-^) + /(a;) -7 - ^TT^^lnA 

93 

+ 24 Lis (1 - a;) - 24 Lis (x) + 12/(a;) + 24C(3) - ^ - IOtt^ 
+ 45S') + C> 



' 2 ^ 



X = 



(19) 



where the function f{x) is given in Appendix. In the formula above besides the logarith- 
mically enhanced terms derived here, we included also the known contribution without 
the large logs, Sq , which is given by Eq. (3) from Ref. [2]. 



4 Numerical Results and Conclusions 



In this way we received the photonic part of the second order next-to-leading logarith- 
mic contribution to Bhabha cross section. The result agree with earlier calculations by 
means of different methods. Our approach allows to get all the next-to-leading contribu- 
tions systematically. It can be applied to any kind of a process, where one has to look 
for the radiative corrections enhanced by large logarithms. In particular, we applied the 
same approach to the description of the contribution of real photon emission to Bhabha 
scattering [to be described elsewhere]. 

Let us compare the numerical values of the leading, next-to-leading, and next-to-next-to- 
leading corrections for two choices of the factorization scale: M = ^/s, which has been 
used in Refs. [1,2], and M = y/^, which has been advocated in Ref. [3]. Since we have the 
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complete answer (19), we can easily choose any other the factorization scale by changing 
the argument of the large logarithm, while the total sum is kept unchanged. For this 
purpose we use the relations 



(1-c) 
t — —s — - — = —xs. 



Ls — Lf 



Inx, 



In 



m 



t 

'2 ■ 



(20) 



In Figures 1 — 5 we show the values of the second order soft and virtual photonic radiative 
corrections in different approximations with respect to the power of the large logarithm. 
Values of the corrections are given in terms of 10^"^ ■ da^°"\ In particular, ^2^] repre- 
sent the leading O (a^L^), the next-to-leading O (a^L^), and the next-to-next-to-leading 
O {a^L^) relative contributions to the cross section, respectively. Since the dependence 
on the parameter A should disappear in the sum of the virtual and soft corrections with 
the remaining three contributions (see Eq. (3)), we put A = 1. In this way we receive 
only an estimate of the magnitude and the relative size of the corrections in different 
approximations. Nevertheless this evaluation helps us to get an idea about the size of the 
unknown second order contributions and to estimate the theoretical uncertainties. 
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Fig. 1. Soft and virtual second order photonic radiative corrections versus the scattering angle 
in degrees for A = 1, ^/s=l GeV; M = ^/s on the left side and M = \/^ on the right side. 

Fig. 1 shows the size of the corrections under consideration relevant for the studies at 
electron-positron colhders of moderate energy (VEPP-2M, BEPS etc.). 

Fig. 2 and Fig. 3 gives us results for the small and large angle Bhabha scattering at 
LEP/SLC, respectively. We checked that for ^/s—200 GeV the plots are very close the 
the ones shown for y^^lOO GeV. 

Figures 4 and 5 give as and idea, how do the corrections behave at higher energies, which 
can be reached at a future linear collider. 



Looking at the plots representing the contributions of different powers in L, we conclude 
that with the proper choice of the factorization scale M = the magnitude of the 

non-logarithmic corrections is below 1 • 10"*^%, everywhere except the region of very large 
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Fig. 2. Soft and virtual second order photonic radiative corrections versus the scattering angle 
in degrees for A = 1, y^=100 GeV; M = ^/s on the left side and M = on the right side. 




Fig. 3. Soft and virtual second order photonic radiative corrections versus the scattering angle 
in degrees for A = 1, ^/s=100 GeV; M = i/s on the left side and M = on the right side. 




Fig. 4. Soft and virtual second order photonic radiative corrections versus the scattering angle 
in degrees for A = 1, ^/s=500 GeV; M = ^/s on the left side and M = on the right side. 

scattering angles {9 > 160°). That region requires a special treatment, and it doesn't 
seem to be of interest for the experiments. Note that the estimate of the size of the non- 
logarithmic second order corrections agree with the one made earlier in Ref. [3]. At the 
same moment it is clear that to reach the 1 ■ 10"*^% level in the precision of theoretical 
description of Bhabha scattering we should take into account the complete O (a^L°) 
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Fig. 5. Soft and virtual second order photonic radiative corrections versus the scattering angle 
in degrees for A = 1, y^=500 GeV; M = ^/s on the left side and M = on the right side. 

calculations including the effects of virtual and real corrections due to pairs and photons. 

Taking the result of the present study and the O (a^L) results of papers [6,7], where the 
pair corrections were evaluated, and of Refs. [16,17], where real photon radiation was taken 
into account, we arrive at the complete result for the second order next-to- leading radiative 
corrections to Bhabha scattering. The results are valid both for the small and large angle 
scattering. To apply the results to data analysis of modern and future experiments at 
electron-positron colliders, we are going to implement them into the Monte Carlo event 
generators LABSMC [18] and SAMBHA [19] for large and small angle scattering, respectively. 
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Appendix. Explicit Formulae for the Corrections 



The differential Born level cross section reads 




(A.l) 



The one-loop virtual and soft corrections read 
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da^ = da^°^'^^{4 In ^ (l - + In (^)) - + 2L, In (i-^) - ln'{x) 

+ ln'(l-x) + 3L,-4 + /(x)}, (A.2) 

+ ln2(x) -ln2(l -x) - — + 2Li2(l -x) -2Li2(a;)j. (A.3) 



/(x) — {1 — X + X 



2\-2 



TT 



^(4 - 8x + 27x^ - 26x^ + 16x^) 



1 1 

+-(-2 + 5x - 7a;^ + 5x^ - 2x^) ln^{l - x) + -x{3 - x - 3x'^ + Ax^) ln^{x) 

+-(6 -8x + Qx"^ - 3x^) \n{x) - -x{l + x'^) ln(l - x) 

2 2 



+-(4 -8x + Ix^ - 2x^) ln(x) ln(l - x) 



(A.4) 



The dilogarithm and the Riemann zeta-function are defined as usually: 



1 

{x) ^ j<i-y 



ln(l — xy) 

y 



OO 1 



n=l 



TV 



(A.5) 
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